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Deep connections are known to exist between scale-free networks and non-Gibbsian statistics. For 
example, typical degree distributions at the thermodynamical limit are of the form P{k) oc , 

where the g-exponential form = [1 -I- (1 — q)z\ i-® optimizes the nonadditive entropy Sq (which, 
for g ^ 1, recovers the Boltzmann-Gibbs entropy). We introduce and study here d-dimensional 
geographically-located networks which grow with preferential attachment involving Euclidean dis¬ 
tances through r~P^ {a a > 0). Revealing the connection with g-statistics, we numerically verify 
(for d =1, 2, 3 and 4) that the g-exponential degree distributions exhibit, for both g and k, universal 
dependences on the ratio OLAjd. Moreover, the g = 1 limit is rapidly achieved by increasing OLAjd 
to infinity. 

PACS numbers: 89.75.He, 05.70.-a, 05.45.Pq, 89.75.Da 


Networks emerge spontaneously in many natural, arti¬ 
ficial and social systems. Their study is potentially im¬ 
portant for physics, biology, economics, social sciences, 
among other areas. For example, many empirical stud¬ 
ies have identified peculiar properties in very different 
networks such as the Internet and online social net¬ 
works (e.g., Facebook), citations networks, neurons net¬ 
works [ll-S|, to quote but a few. An ubiquitous class of 
such networks is constituted by the scale-free ones (more 
precisely, asymptotically scale-free). As we shall soon 
verify, these networks can be seen as a particular ap¬ 
plication of nonextensive statistical mechanics, based on 
the nonadditive entropy Sq = {q G TZ; Si = 

Sbg= — Pi In Pi, where BG stands for Boltzmann- 
Gibbs) [3,0- This current generalization of the BG en¬ 
tropy and corresponding statistical mechanics has been 
widely successful in clarifying the foundations of thermal 
statistics as well as for applications in complex systems 
in high-energy collisions at LHC/CERN (GMS, ALICE, 
ATLAS and LHCb detectors) and at RHIC/Brookhaven 
(PHENIX detector) 0 , cold atoms [ 3 ], dusty plasmas 0 , 
spin-glasses 0, trapped ions 0, astrophysical plasma 
[Tl| , biological systems fT3. t voe-II superconductors [10 , 
granular matter [10 (see [l5j|). 

The deep relationship between scale-free networks and 
g-statistics started being explored in 2005 [TMl , and 
is presently very active [l9l - l^ . The basic connec¬ 
tion comes (along the lines of the BG canonical en¬ 
semble) from the fact that, if we optimize the func¬ 
tional Sq[P{k)] = .^^ith the constraint 

(k) = f dkkP(k) = constant or analogous {k being 
the degree of a generic site, i.e., the number of links 
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that arrive to a given site; P{k) denotes the degree or 
connectivity distribution), we straightforwardly obtain 
P{k) = P{0)eq'^^'^ = P(0)/[1 + {q — l)fc/K]^, which 
turns out to be the generic degree distribution for virtu¬ 
ally all kinds of scale-free networks. The g-exponential 
function is defined as = [1 -I- (1 — g)z]i^ (ef = e^). 
We verify that, for g > 1 and k oo, P(^ ~ l/fc'*' with 
7 = l/(g—1). The classical result 7 = 3 [^ corresponds 
to g = 4/3. 

In the present work we address the question of how 
universal such results might be, and more specifically, 
how P{k) varies with the dimension d of the system? 

Our growing model starts with one site at the origin. 
We then stochastically locate a second site (and then 
a third, a fourth, and so on up to N) through the d- 
dimensional isotropic distribution 

P(d oc (ac > 0; d= 1,2,3,4), (1) 

where r > 1 is the Euclidean distance from the newly 
arrived site to the center of mass of the pre-existing sys¬ 
tem (in one dimension, r = \x\; in two dimensions, r = 
-\- y'^; in three dimensions r = x^ + -\- and 

so on); p{r) is zero for 0 < r < 1; the subindex G stands 
for growth. We consider ac > 0 so that the distribu¬ 
tion P{r) is normalizable; indeed, ^d-ij,-(d+aG) = 

dr l/r^+“'^, which is finite for ao > 0, and diverges 
otherwise. See FigH] 

Every new site which arrives is then attached to one 
and only one site of the pre-existing cluster. The choice 
of the site to be linked with is done through the following 
preferential attachment probability: 

r. j,. .—OA 

= €[0,1] (aA>0), (2) 

where ki is the connectivity of the i-th pre-existing site 
(i.e., the number of sites that are already attached to site 
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i), and Tij is the Euclidean distance from site i to the 
newly arrived site j ; subindex A stands for attachment. 

For a A approaching zero and arbitrary d, the physical 
distances gradually loose relevance and, at the limit a a = 
0 , all distances becomes irrelevant in what concerns the 
connectivity distribution, and we therefore recover the 
Barabasi-Albert (BA) model [IJ , which has topology but 
no metrics. 
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FIG. 1. Distribution of A = 500 sites obtained with Eq. © 
for aA = 2.0, ac ~ 0.0, and d = 1, 2, 3. 
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FIG. 2. Connectivity distribution for d = 1,2, 3, 4, aA ~ 2.0 
and typical values for ac- The simulations have been run 
for 10^ samples of A = 10® sites each. We verify that P{k) 
independs from aa (Vd). 

Large-scale simulations have been performed for the 


(d = 1 , 2 ,3,4) models for fixed {ac,o:A), and we have 
verified in all cases that the degree distribution P{k) is 
completely independent from ac- see Fig. O Using this 
fact, we have arbitrarily fixed ac = 2 , and have nu¬ 
merically studied the influence of (d, a^) on P{k): see 
Figs. [3] and © In all cases, the ^-exponential fittings 
P{k) = P(0)eg with q > 1 and k > 0 have been 
remarkably good. The best htting values for {q, k) are 
indicated in Fig. [S] From normalization of P(k), P{0) 
can be expressed as a straightforward function of ( 9 , k). 

Our most remarkable results are presented in Fig. [51 
namely the fact that both the index q and the character¬ 
istic degree (or “effective temperature”) k do not depend 
from {aA,d) in an independent manner but only from 
the ratio aAfd. This nontrivial fact puts the growing 
d-dimensional geographically located models that have 
been introduced here for scale-free networks, on simi¬ 
lar footing as long-range-interacting many-body classi¬ 
cal Hamiltonian systems such as the inertial XY planar 
rotators [ 2 ^ (possibly the generic inertial n-vector rota¬ 
tors as well [ 2 ^) and Fermi-Pasta-Ulam oscillators, 
assuming that the strength of the two-body interaction 
decreases with distance as l/(dfstance)“. Moreover, as 
first pointed out generically by Gibbs himself [ 2 ^, we 
have the facts that the BG canonical partition function 
of these classical systems anomalously diverges with size 
for 0 < a/d < 1 (long-range interactions, e.g., gravita¬ 
tional and dipole-monopole interactions) and converges 
for a/d > 1 (short-range interactions, e.g., Lennard- 
Jones interaction), and the internal energy per particle is, 
in the thermodynamical limit, constant for short-range 
interactions whereas it diverges like for long- 

range interactions, N being the total number of particles. 

If all these meaningful scalings are put together, we 
obtain a highly plausible scenario for the respective do¬ 
mains of validity of the Boltzmann-Gibbs (additive) en¬ 
tropy and associated statistical mechanics, and that of 
the nonadditive entropies Sq (with q ^ 1) and associated 
statistical mechanics. 

Finally, we notice in Fig. ©that both q and k approach 
quickly their BG limits {q = 1) for a A/d 00 . More¬ 
over, the same exponential appears in both heuris¬ 

tic expressions for q and k. Gonsequently, the following 
linear relation can be straightforwardly established: 

K ~ 4.90- 3.45 g. (3) 

In fact, this simple relation is numerically quite well sat¬ 
isfied as can be seen in Fig. [71 Its existence reveals an 
interesting peculiarity of the nature of g-statistics. If in 
the celebrated BG factor ^-energy/kT^ corresponding to 
g = I, we are free to consider an arbitrary value for T, 
how come in the present problem, k is not a free param¬ 
eter but has instead a fixed value for each specific model 
that we are focusing on? This is precisely what occurs in 
the high-energy applications of g-statistics, e.g., in quark- 
gluon soup where g = 1.114 and T = 135.2 Meu, as 
well as in all the LHG/GERN and RHIC/Brookhaven 
experiments Q. Another example which is reminiscent 
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FIG. 3. Degree distribution for d = 1 {blue diamonds), 
2 {green triangles){magenta squares), 4 {grey circles), 
and typical values of a a, with aa = 2.0. The simnlations 
have been run for 10® samples of = 10® sites each. 

of this type of behavior is the sensitivity to the initial 
conditions at the edge of chaos (Feigenbaum point) of 
the logistic map; indeed, the inverse ^-generalized Lya- 
punov exponent satisfies the linear relation 1/Ag = 1 — q 
|30l| . The cause of this interesting and ubiquitous feature 
comes from the fact that g-statistics typically emerges 
at critical-like regimes and is deeply related to an hier¬ 
archical occupation of phase space (or Hilbert space or 
Fock space), which in turn points towards asymptotic 
power-laws (see also [3l|). In other words, k plays a role 
analogous to a critical temperature, which is of course 
not a free parameter but is instead hxed by the specific 
model. 
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FIG. 4. Fittings of the d = 1,2,3, 4 connectivity distribu¬ 
tions with the function P{k) = P{0)eq~'^^*° , where e, = 
[1 + (1 — The data are those of Fig. [S] Top: 

log-log representation. Bottom: lnq[P{k) / P{Q)] versus k rep¬ 
resentation. The fitting parameters are exhibited in Fig. [5] 
The numerical failure, at large enough values of k, with regard 
to straight lines are finite-size effects that gradually disappear 
when we approach the thermodynamic limit N —>■ oo. 
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FIG. 5. q and n for d = 1,2,3, 4. For oa = 0 and Vd, we 
recover the Barabasi-Albert universality class q = 4/3 (hence 
7 = 3) [2J, which has no metrics. 
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